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Abstract 

After a brief introduction we review the nonperturbative weak noise 
approach to the KPZ equation in one dimension. We argue that the strong 
coupling aspects of the KPZ equation are related to the existence of local- 
ized propagating domain walls or solitons representing the growth modes; 
the statistical weight of the excitations is governed by a dynamical action 
representing the generalization of the Boltzmann factor to kinetics. This 
picture is not limited to one dimension. We thus attempt a generalization 
to higher dimensions where the strong coupling aspects presumably are 
associated with a cellular network of domain walls. Based on this picture 
we present arguments for the Wolf-Kertez expression z — (2d+l)/(d+l) 
for the dynamical exponent. 
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1 Introduction 

Driven systems far from equilibrium constitute an enormous class of natural phe- 
nomena including turbulence in fluids, interface and growth problems, chemical 
reactions, processes in glasses and amorphous systems, biological processes, and 
even aspects of economical and sociological structures. 

In recent years much of the focus of modern statistical physics and soft con- 
densed matter has shifted towards such systems. Drawing on the case of static 



'permanent address 



and dynamic critical phenomena in and close to equilibrium, where scaling and 
the concept of universality have successfully served to organize our understand- 
ing and to provide a variety of calculational tools, a similar strategy has been 
advanced towards the much larger class of nonequilibrium phenomena with the 
purpose of elucidating scaling properties and more generally the morphology or 
pattern formation in a driven nonequilibrium state. 

In this context the Kardar-Parisi-Zhang (KPZ) equation, describing the 
nonequilibrium growth of a noise-driven interface, provides a simple contin- 
uum model of an open driven nonlinear system exhibiting scaling and pattern 
formation The KPZ equation for the time evolution of the height h{x,t) 
has the form |^ 

dh(x t) A -< -< 

— = i/V2/i(f , t) + -Vh[x, t) ■ Vh[x, t) + T]{x, F , (1) 
{tj{x, t)r]{x', t')) = A(5(f - ^)5{t - t') . (2) 

Here v characterizes the damping term, A controls the growth term, is a gen- 
eral drift term, and ri[x,t) a locally correlated white Gaussian noise modelling 
the stochastic nature of the drive or environment; the noise correlations are 
characterized by the noise strength A. 

Starting from a given initial interface the noise drives after a transient period 
the interface into a stationary fluctuating state. Choosing a co-moving frame by 
setting F — (A/2)(V/i • Vh) the mean height (h) decays to zero but the fluctua- 
tion spectrum about (h) = 0, characterized by the moments (/i") and generally 
by the distribution P({h{x)}, t), assumes a nontrivial and largely unknown form. 
In terms of the local slope u = Vh the KPZ equation assumes the form of the 
Burgers equation driven by conserved noise (^d/dt — [u ■ V)^ u = vV^u + V77 
which within the framework of perturbative dynamic renormalization group 
(DRG) theory was first studied as a model of noise-driven irrotational fluid flow 
. This analysis was reviewed and extended considerably in the context of the 
KPZ equation [||. The DRG analysis basically addresses the long time - large 
distance scaling properties as illustrated for example by the dynamic scaling 
hypothesis [Q applied to the stationary height-height correlations 

{{h{x,t)-h{x',t')f) = \x~x'\^^G{\t-t'\/\x~x'\') . (3) 

Here the roughness exponent C, the dynamic exponent z, and the scaling func- 
tion G together define the KPZ universality class. Power counting or a DRG 
analysis identifies c? = 2 as the lower critical dimension. An e-expansion about 
d = 2 yields a kinetic phase transition at a finite growth strength A above d = 2 
from a weak coupling smooth phase with z = 2 and C = to a strong coupling 
rough phase with largely unknown scaling exponents. 

In d = 1 the DRG analysis fortuitously yields z = 3/2 and ( = 1/2. This 
result also follows from the Galilean invariance of (^, x — > x — Xuot, u — > u + uq, 
/i — > h + uq ■ X, which since A is unrenormalized under the DRG, implies the 
scaling law z + ( = 2. Together with the known stationary distribution, Pg^ = 
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exp(— (z//A) J dx u^) Q, of independent slope fluctuations, yielding C, ~ 1/2, 
we then obtain z = 3/2. However, unlike critical phenomena, the DRG analysis 
even carried to orders beyond first loop order [Q fails to yield detailed insight 
into the scaling and pattern formation mechanism in the KPZ equation. 

A mode coupling analysis (MC) of the KPZ equation has also been carried 
out both with regard to a determination of the scaling function G vn d — 1 
and in search for an elusive upper critical dimension assumed to be at d = 4 
beyond which z = 2 and C = Hi ■ Of other approaches to the scaling properties 
associated with the KPZ universality class we mention briefly i) the mapping of 
the KPZ equation to the problem of directed polymers in a quenched random 
environment Q , a problem in ill-condensed matter, ii) exact results obtained for 
the d = 1 asymmetric exclusion model which is in the same universality class 
as the KPZ equation, and finally iii) recent exact results for a discrete d = 1 
polynuclear growth model in the KPZ universality class through the association 
with random matrix theory 0|. In Fig. |l| we have in a plot of the renormalized 
coupling strength versus the dimension summarized the scaling properties of the 
KPZ equation. 

In summary, the KPZ equation as a continuum model of an intrinsic nonequi- 
librium problem and its relationship to problems in discrete driven lattice gases, 
random systems of the spin glass type, and random matrix theory, is of broad 
and paradigmatic interest. 



2 Strong coupling theory 

In recent work we have proposed yet another approach to the scaling and pattern 
formation of the KPZ equation |^ which we review below. The strong coupling 
approach is developed by focusing on the Fokker-Planck equation 

^ dP{h{x),t) ^ ^ H = J dx n{x) , (4) 

associated with the KPZ equation (^. Here the Hamiltonian or Liouvillian H 
drives the transition probability P{h{x, t)) for the height profile h{x) at time t. 
The Fokker-Planck equation has the form of a Schrodinger equation in imaginary 
time for the real positive "wave function" P{h{x),t). The noise strength A 
plays the role of an effective "Planck constant". Since the KPZ equation is 
a stochastic field equation the associated Fokker-Planck equation is a vastly 
more complicated functional-differential equation. However, applying the WKB 
approximation valid in the "semiclassical" weak noise limit A ^ by setting 
(hi in the initial height configuration) 

S{hi{x) ^ h(x),t)] 



P[hi{x) — > h{x),t)] (X exp 
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Figure 1: We summarize the scaling properties of the KPZ equation in a plot 
of the renormalized coupling constant (the fixed point value) AA^/i/^ versus the 
dimension d of the system. Below d — 2 the scaling properties are determined 
by the strong coupling fixed. Above d — 2 the system exhibits a kinetic phase 
transition from a weak coupling phase to a strong coupling phase. We have 
also indicated three conjectures for the exponent z: The operator expansion 
conjecture zl by Lassig Q and the conjectures based on numerics, zwk by 
Wolf and Kertesz and zkk by Kim and Kosterlitz . 
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a "principle of least action" becomes applicable and the "weight function" or 
action takes the form, 



fh{x),t 

S{hi{x) h{x),t) = j dxdt 



PixMf-^-nx) 



'hi{x),t=0 

This implies the "classical" Hamilton equations of motion 



(6) 



^IM^ ^ y\j^h{x, t) + ^Vh[x, t) ■ Vh{x, t)-F + p{x, t) , (7) 
^ -j/V2p(f , t) + AVp(f , t) • V/t(f , t) + \p(x, t)V^h(x, t) , (8) 

at 

derived from the Hamiltonian density 

vV'^h{x, t) + ^Vh{x, t) ■ Vh{x, t) - F + ^p{x, t) 



TC — p{x, t) 



(9) 



The field equation and (||) are of the difi^usive-advective type and replace 
the KPZ equation. The deterministic "noise field" p{x, t) corresponds to the 
stochastic noise ri{x,t) in (||). 

The prescription for the analysis of the KPZ equation is now straightforward. 
In order to determine the transition probability from an initial height profile 
hi{x) at time t = to a final height profile h{x) at time t we must evaluate 
the "classical" action associated with a "classical" orbit from hi{x) to h{x) 
traversed in time t. The strong coupling features enters in the determination 
of the relevant orbits which should be nonperturbative solutions of the field 
equations. As is well-known from the WKB approximation in quantum theory 
and quantum field theory, the weak noise limit has no direct bearing on the 
coupling strength. The action in the present nonequilibrium context plays the 
same role as the energy in the Boltzmann factor in equilibrium. 

Since the "classical" system is conserved the orbits are confined to lie on 
energy surfaces H ~ const. Here the zero energy surface H = plays a central 
role in determining the stationary distribution Pg^f^ih) — limt^oa P{hi —^ h,t). 
Unlike the case in mechanical systems, the zero energy surface is unbounded and 
has a characteristic submanifold structure reflecting the underlying stochastic 
nature of the KPZ equation. The transient submanifold for p — corresponds 
to the damped noiseless KPZ equation for ?/ = 0, whereas the in general un- 
known stationary submanifolds is "orthogonal" to p = in the sense that the 
J dxTi = with Ti. given by (^) . The submanifolds intersect in the saddle point 
{h,p) = (0, 0) corresponding to ergodic behavior. At long times a specific orbit 
is unbounded; however, an orbit from hi to h in time t is finite. At long times 
the orbit migrates from the finite energy surface to the zero energy submani- 
folds. After an initial transient period the noise drives the orbit close to the 
saddle point. Ergodic behavior is established and the orbit moves out along 
the stationary submanifold. At long times the orbit converges to the stationary 
submanifold and the stationary stochastic state is attained. This generic phase 
space behavior is depicted schematically in Fig. |^. 
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Figure 2: We depict the generic phase space behavior of the WKB approxima- 
tion apphed to the Fokker-Planck equation. The stationary saddle point (SP) 
is at the origin. The transient zero energy submanifold and the stationary zero 
energy submanifold are indicated. At long times the orbit from hi to h in time 
t converges to the submanifolds and passes close to the saddle point. 



3 Morphology and scaling in one dimension 

The strong coupling WKB analysis becomes relatively simple in one dimension 
[pi . The stationary zero energy submanifold is given by p{x,t) = -2vV'^h{x,t), 
yielding the stationary distribution [Q Pf^i{h) oc exp[— (i//A) / dx{dh/dx)'^]. 
Moreover, the field equations (0) and (^) admit simple cusp solutions, in the 
static case of the form h{x) — {2v/X) log | cosh{{X/2i/)u{x — xo)) \ centered about 
Xq and for large x approaching the constant slope form h ^ u\x — xo\. In terms 
of the local slope dh/dx the cusps correspond to solitons or domain walls of the 
kink-like form dh/dx — utanh((A/2:/)|u|(a; — xq)). The right hand soliton for 
u > moves on the p — submanifold and is the well-known viscosity-smoothed 
shock wave solution of the noiseless Burgers equation {d/dt — AuV)it = vV'^u. 
The left hand soliton for u < moves on the p = —2v'S/'^h submanifold and 
is a solution of the growing Burgers equation {d/dt — AuV)m — —vV^u with 
negative damping constant. The cusps or solitons form the elementary exci- 
tations or growth modes in a heuristic many body description of a growing 
interface. Boosting the cusps or solitons to a finite propagation velocity v by 
means of the Galilean transformation above we obtain the kinematic condition 
V = — (2/A)(m+ + U-), where u+ and U- are the right and left boundary values 
of the solitons, respectively. Matching a set of right and left hand solitons ac- 
cording to the kinematic condition we thus obtain a many body or multi-soliton 
representation of a growing interface. 

Imposing periodic boundary conditions it is easily seen by inspection that 
the propagation of solitons lead to a growing interface. A particular simple 



6 



h 




ih 


ih 






DOWN 
CUSP 


UP 
CUSP 




m 






dh 

dx 


dh 


dh 

RHS 1 


LHS^ 









Figure 3: We depict the single cusp and soliton modes and a multi soliton 
representation of a growing interface. We show both the height profile and the 
slope profile. 

growth mode corresponds to the propagation of a pair of matched right and left 
hand solitons. This mode is equivalent to the addition of layer upon layer in 
the evolution of the height profile. In Fig. ^we have shown the single cusps and 
solitons and a multi soliton representation of a growing interface 

The WKB approach also allows us to associate a dynamics with the propa- 
gation of growth modes. Only the noise- induced left hand solitons are endowed 
with dynamical properties. In terms of the boundary values u+ and U- we ob- 
tain the Galilean invariant action S = (l/6)Ai'|u+ — from the Hamiltonian 
H the energy E = (2/3)Ai^(ui]_ — v?_\ and from 11 — J dxhS/p the momentum 
U — ^[u^—u^). Considering specifically a two-soliton mode composed of a right 
hand and left hand soliton with amplitude u and with vanishing boundary con- 
ditions it carries action 5* = (l/6)|upt, energy E — — (2/3)Ai^|up, momentum 
n = and propagates with velocity v — —Am/2. Eliminating the velocity 

we obtain the soliton dispersion law and from (|^) the pair soliton distribution, 
X denotes the center of mass coordinate, 

E = -l^\^\'^' and F(a:,t)cxexp(-1^^) . (10) 

From these two expressions we draw two results: i) It follows easily from the 
spectral properties of the height correlations that the exponent 3/2 in the soliton 
pair dispersion law can be identified with the dynamic exponents z. The dynam- 
ical scaling is thus associated with the low lying gapless soliton modes, ii) In 
the stochastic growth mode characterized by a gas of independent pair solitons 
the solitons perform random walk with the characteristics of superdiffusion, e.g., 
the distribution implies the mean square displacement {x^){t) cx (AA^/j^)^/^t^/^ 
with z = 3/2. 
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4 Towards higher dimensions 



In dimensions beyond d = 1 the scaling and morphology are more involved. 
From the DRG and the mapping to directed polymers if follows that d = 2 
is a lower critical dimension. Above d = 2 the system undergoes a kinetic 
phase transition at a finite coupling from a weak coupling smooth phase with 
exponents z = 2 and = to a strong coupling rough phase. An operator 
expansion method yields z — Ad/{2d+ 1) for d > 1 Two other heuristics 
proposals are z = 2(d + 2)/(d + 3) and z = {2d + l)/{d + 1) based on 
numerics; they both yield z = 3/2 in cZ = 1. Common to all these expressions is 
that z ^ 2 for d — > oo, i.e., yielding an infinite upper critical dimension. 

Within the present WKB approach the evolution of the morphology of the 
growing interface is given by the appropriate strong coupling solutions of the 
"classical" equations and (||) . In one dimension we were able analytically to 
find soliton solutions and thus "parameterize" a growing interface in terms of 
a gas of kinematically matched solitons. We have also analyzed the equations 
numerically in one dimension but here the negative diffusion coefficient in the 
equation for the noise field p renders the coupled equations numerically unstable 
and we were only able to verify certain soliton collision configurations. The same 
instability problem will remain in higher dimension and it is at the moment not 
clear how to find solutions numerically. 

Starting from random initial conditions the time evolution of the noiseless 
Burgers equation is damped [ p2[ . Since the nonlinear Cole-Hopf transformation 
[0 relates the slope field to a linear diffusive field the transient morphology can 
be analyzed by saddle point methods and is characterized by a gas of localized 
right hand solitons connected by ramp solutions; for the height field this pat- 
tern formation corresponds to a gas of downward cusps connected by parabolic 
segments. In two and higher dimensions one finds a corresponding morphology 
[ |l2| . The height field forms a damped cellular network of growing "domes" 
connected by "valley's" . The domes and valleys correspond to the parabolic 
segments and downward cusps in one dimension, respectively. 

When the noise is "turned on" the right hand soliton in one dimension is sup- 
plemented by a noise induced left hand soliton in the WKB interpretation and 
as discussed above the morphology can be discussed in terms of connected right 
hand and left hand solitons. For the height field one version of this morphology 
corresponds to growing and decaying plateaus. The moving steps connecting 
the plateaus correspond to soliton pairs in the slope field. Similarly, we be- 
lieve that a corresponding morphology takes place in higher dimensions. In for 
example two dimensions the stationary growth morphology in the height field 
is composed of growing and decaying islands forming a cellular network. The 
islands are connected by ramps or steps propagating with a velocity dictated by 
the Galilean invariance of the corresponding Burgers equation. The statistical 
weight of a morphology is then given by the action associated with a given time 
evolution from hi to in a time t. If we parameterize the islands by linear 
segments and ignore vertex corrections we can make use of the one dimensional 
soliton solutions and connect the island by means of pair-solitons of amplitude 
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Figure 4: We depict the cellular network of growing and decaying islands in the 
height profile. The network is formed of domain walls or solitons constituting 
the fundamental growth modes. 



u. The steps will then propagate with a velocity v or order Xu and carry an 
action of order Xvu^t per. unit length. In Fig. ^ we have depicted a growth 
configuration for the height field in two dimensions. 

We now have enough elements of a future more detailed theory to venture a 
proposal for the dynamic exponent z. First we notice that the localized growth 
modes lie on the stationary manifold H — 0. In one dimension this manifold 
is given hy p = —2i'\7^h, yielding the action S oc Ai^jupt and the stationary 
distribution Pgiih) oc exp[— (i^/A) / dx{dh/ dx)"^]. In higher dimensions the sta- 
tionary manifold is unknown and the stationary distribution is expected to be 
non-Gaussian and skew Q. However, assuming that the Galilee invariant ac- 
tion for a step in /i is a function of the slope change 5u (ignoring the vector 
character of u) and expressing S in the form S cx Xv\u\'^P^~^t, where the factor 
arises from the extended character of the higher dimensional "solitons" 
or "domain walls" providing the growth modes, a natural choice interpolating 
between d = 1 and d = oo is to set a = d + 2. Finally, noting that 5u oc v 
and balancing the space and time dependencies of v and I we infer a dynamic 
exponent in agreement with Wolf and Kertez [jll] 

This expression agrees with z = 3/2inc?=l. Inff = 2 we obtain z = 5/3 which 
is close to the numerical value For d ^ cx) we have z — > 2 and the analysis 
yields an infinite upper critical dimension. 
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